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Abstract 

The large spacing phase of the inhnite random matrix chain, which represents the strongly 
conpled two-dimensional 0(2) model on a random planar lattice, is explored. A class of 
solntions valid for large lattice spacings is constrncted. It is proved that these solntions 
exhibit the critical exponents characteristic of pnre two-dimensional gravity. The char¬ 
acter expansion for the chain model is developed and an order parameter governing the 
Kosterlitz-Thonless phase transition is identihed. 

1. Introduction 

String theories with the matter central charge one are of special interest. They are 
positioned in between the exactly solved c < 1 minimal models interacting with qnantnm 
gravity |]l| and the yet mysterions c > 1 noncritical string theories 0 0 . 
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The c = 1 string where the target space is an inhnite one-dimensional line can also be 
analyzed exactly 0. However, new problems arise once the target space is compactihed. In 
that case the theory has been solved only for snfhciently large values of the compactihcation 
radius R > Rcr- At R = Rcr the c — 1 model undergoes a phase transition, and there is 
evidence that for smaller R it behaves like a theory of pure two-dimensional gravity, without 
any matter at all 0 [^ . 

The transition at R = Rcr is quite similar in nature to the famous Kosterlitz-Thouless 
phase transition in the two-dimensional 0(2) model 0. Naively, these transitions are not 
automatically present in continuum theory and appear only if the worldsheet is discretized. 
Both of them are induced by topologically nontrivial held conhgurations—vortices—where 
the compactihed string, or a two-component unit vector in the 0(2) model winds around a 
unit circle as we follow the boundary of an elementary worldsheet plaquette. The dynamical 
role of vortices in the theory depends on how the vortex energy and entropy compare to 
each other. As a result, while the ehect of vortices is negligible for R > Rcr (or in the 
weakly coupled phase of the 0(2) model) the vortices actually dominate the dynamics at 
R < Rcr (respectively, in the strong coupling phase of the 0(2) model.) Furthermore, the 
well known R a'/R duality of one-dimensional compactihed string theory holds only if 
the vortices are completely ignored. It is possible to show |P that the contributions to the 
free energy coming from the vortexless sector of the c = 1 theory do respect the T-duality, 
whereas the vortex contributions manifestly break it. 

To understand the dynamics of vortices interacting with quantum gravity is a long¬ 
standing problem . Technically, this problem can be formulated in a very natural way 
using the language of random matrix models. Indeed, there the two-dimensional world- 
sheet is discretized by construction. Thus any matrix model describing the compactihed 
one-dimensional string theory will contain vortices from the very beginning. The difficulty 
is, all such models involve an inhnite number of interacting matrices and, because of that, 
have not yet been solved. 

In this paper we shall explore the simplest among these models, a one-dimensional 
inhnite random matrix chain dehned by the partition function H H 


1C 


Z = 


/ -|-oo I' -hCXD 

n dMr, exp|-iVTr ^ 

n — — ri= — oo - 


{Mn+l - 
2e 


+ eH(M^) 
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Here all are N x N Hermitian matrices while V{M) is a polynomial potential such as, 
for example, m^M^/2 -|- gM^/3 or m^M^/2 -|- gM'^/4. 

It is possible to show |]^ that the large N limit of the chain model does indeed describe 
the one-dimensional bosonic string theory compactihed on a circle of radius R = 1/e or. 








equivalently, the two-dimensional 0 (2) nonlinear sigma model coupled to quantum gravity. 
That is to say, the leading term in the large N expansion of Z summarizes correctly the 
vortex properties on topologically spherical random surfaces. At the same time the 1 /N 
corrections to the chain partition function do not represent the higher genus compactihed 
c = 1 amplitudes. Instead, they correspond to a different theory—a one-dimensional 
string with discretized target space. On a spherical worldsheet this theory is related to the 
compactihed one-dimensional string by a duality transformation, and the string partition 
functions for the two coincide. 

It has long been known that the inhnite matrix chain undergoes a Kosterlitz- 
Thouless phase transition at a certain e = Ccr- Throughout the region e < Ccr where 


vortices can be neglected the critical indices of (|1.1| ) do correspond to c = 1. On the other 
hand, very little, if anything, is known about the properties of the matrix chain for e > Ccr. 
Quite remarkably, it is this regime, totally dominated by vortices, that is the hardest to 
analyze by traditional means. 

Certainly, in the limit of inhnitely large e the matrix chain partition function decouples 
into a product of independent one-matrix models. Therefore, one might conjecture that, 
at least for large e, the critical indices of Z are those of pure two-dimensional gravity. We 
shall see that such a conclusion is in fact true for all e > Ccr- As a consequence, the effects of 
dimensional reduction in string theories and held theories may in fact be much more sim¬ 
ilar than one usually believes. Consider, say, a held theory on a manifold with a compact 
dimension. When the compactihcation radius vanishes the compactihed dimension ehec- 
tively disappears, and we obtain a (perhaps modihed) held theory in one less dimension. 
In string theory the picture is quite diherent. There because of the R •«-*• a'/R duality the 
small and large compactihcation radii appear to be strictly equivalent, and taking R to 
zero does not reduce the ehective dimension. However, once the nonperturbative ehects— 
such as vortices—are included, the situation may change drastically and the parallels with 
held theory may get restored. Exactly that happens in our one-dimensional example where 
for R < Rcr the compactihed string theory behaves as if the one-dimensional target space 
circle had completely disappeared and the target space dimension was equal to zero. 

Below we shall construct, for certain V (M), an exact solution of the matrix chain valid 
at e > Ccr- Remarkably, the computational tools needed for that are quite interesting on 
their own. As we shall see, the large N asymptotics of the matrix chain model is related to 
a certain one-dimensional hydrodynamic system. The Kosterlitz-Thouless phase transition 
corresponds then to the formation of a shock-type conhguration in the moving huid. 



To state this more precisely, let us first rescale all matrices = 
disappears from the kinetic term in the partition function, 

/ oo ^ oo 
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the new potential U{A4) being related to V{M) by 
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V(M) = + tV{^€M). 


(1.3) 


PiX N ^ OO the eigenvalues Ai^ri ■ ■ ■ Aiv,n of any given matrix M.n in (|1.2|) condense to form 
a smooth distribution Pn(A). Actually, due to the translational invariance of the matrix 
chain this distribution is the same for any chain site, Pn(A) = p(A). It turns out that the 
density p(A) can be found as a solution to the following hydrodynamic problem []TT| [p!2| . 
Consider a one-dimensional droplet of compressible fluid with a special equation of state 
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P = - p- 

3 ^ 


(1.4) 


where P is the local pressure and p —the local fluid density. Imagine that at time t = 0 
the density p equals (the unknown) p{x) and the initial fluid velocity is 


v{x) = -U'{x) — X. 


(1.5) 


Now, demand that after one unit of time, at t = 1, the density of the droplet returns to 
its initial value, 

p(x, t = 1) = p(x, t = 0) = p{x). (1-6) 


This condition hxes p{x) once U{x) (and hence n(a:)) is given, so that the problem of the 
matrix chain reduces to solving the hydrodynamic Euler equations. 

As a consequence of ( p..6|) , together with the time reversal symmetry of Euler equa¬ 
tions, the initial and hnal fluid velocities are opposite, v{x,t = 1) = —v{x^t = 0) and, 
furthermore, v[x^t = 1/2) = 0. The density at t = 1/2, which we denote pi/2(a:), shall 
play a very important role in our analysis. Its value at the origin = pi/2(0) provides an 
“order parameter” for the Kosterlitz-Thouless transition. That is to say, C = 0 ^my 
e < Ccr but C > 0 everywhere in the phase of large e. 

Of course, the Kosterlitz-Thouless transition is not explicitly associated with the 
spontaneous breaking of any symmetry and using the term “order parameter” may appear 
purely formal. However, this phase transition does reflect a nonanalytic change in the 
group representation structure of (|1 . 1|) at large N. In section 4 we shall demonstrate that 
the integrand of Z can be expanded as a series over characters of SU{N). At inhnite N only 
one representation in this series, with the Young tableau row lengths (ni,..., njv) all of 
order N, makes an important contribution. Such representations are usually characterized 
by the Young tableau density pi{h) dehned as the density of points hi = (rii — i)/N -|- 1 in 
a small interval around h. Since the row lenghts are all ordered, pi{h) can never be greater 
than one. 






Fig. 1 : The Young tableaus for the representations dominant at e < 1 (left) and 
e > 1 (right). 


Whenever the Young tableau density saturates this upper bound, pi{h^) = 1, the 
corresponding Young tableau develops a gap. It turns out that such gap formation is just 
another aspect of the same Kosterlitz-Thouless phase transition at e = Ccr- Moreover, the 
“order parameter” ( = pi/ 2 ( 0 ) happens to be very directly related to the gap size h*, 

C = Pi/2(0) = -V^. (1.7) 

TT 


Such a correspondence unihes, at least technically, the Kosterlitz-Thouless transition and 
the Douglas-Kazakov phase transitions known from two-dimensional QCD and the dually 
weighted graph models [^|^. The analogy between those transitions goes, in fact, rather 
far. Dynamically, all of them are induced by the topologically nontrivial classical solutions 
of underlying theories—vortices in the 0(2) model and instantons in QCD 2 respectively, 
which become statistically dominant at strong coupling |jT^ . 


Below we shall hrst explore the large e limit of the chain model ( p..l|) . This shall 
be done in the next section. In section 3 we shall construct a family of exact solutions 
describing the matrix chain in the interval e > Ccr- The character expansion methods shall 
be presented in section 4. Finally, in section 5 we shall conhrm our conclusions by studying 
an explicit 1/e expansion for the case of a pure quartic potential U{x). 




























2. The Infinite Matrix Chain at Large Lattice Spacings 


To begin with, we shall demonstrate how the hydrodynamic representation (|1.4| )- 
( p..6|) reprodnces the solntion of the infinite matrix chain in the limit e 3> Ccr- We shall 
nse this special case to understand certain properties of the model—especially the analytic 
structure of the midway eigenvalue density Pi/2{x) —which turn out to hold more generally 
for any finite e > Ccr- 

For large e the “kinetic” term in the matrix model action, tr(M„_|_i — M„)^/2e is 
much smaller than the potential energy term etrld(Mn). Equivalently, one can neglect the 
nearest neighbor coupling term tr(AfnWln+i) in (|1.2|) thereby reducing Z to a product of 
identical one-matrix partition functions 


Z = 
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( 2 . 1 ) 


The eigenvalue density p(A) for the factorized ensemble (p.l|) is the same as the density in 
a one-matrix model with the potential U (Af) and can be found from the Riemann-Hilbert 


equation 16 


f,(X')dX' 


( 2 . 2 ) 


For concreteness let us consider the specific case of a quartic double-well potential 


2 ~ 

V{M) = (2.3) 

with < 0 and ^ > 0. By rescaling M, g and the lattice spacing e it is always possible to 
make the absolute value of anything we want. From now on we shall choose — —4. 
As we shall see later, for this the position of the Kosterlitz-Thouless point is fixed at 
Ccr = 1- The potential U corresponding to such V has the form 


U{x) 


2 2 4 

11 X gx 

2 ^ ~ 


(2.4) 


with = 2(2e^ - 1) and g = ge^. 

Strictly speaking, the one-matrix model with a double-well V{M) describes not the 
c = 0 but c = —2 quantum gravity. To represent the pure two-dimensional quantum 
gravity one must change the overall sign of V taking m? > 0 and g < 0. The computational 
methods presented below are equally applicable in both cases but the formulas happen to 
be simple for the double-well model. Therefore we shall discuss the double-well c = —2 
theory first and then outline the modifications arising when c = 0. 









For the potential U{M) given by (2^) the solntion of the Riemann-Hilbert equation 
I) is quite easy to hnd directly. At sufficiently large g (that is, in the strong coupling 
phase) such solution is given by 

7 rp{x) = {Ax‘^ + B)^/l — (2.5) 


where the coefficients A, B and C equal 

+ 12g - 

B=-^{Vg^ + 12g-2g^) 

A = ^ 

2C' 


( 2 . 6 ) 


Equation (|2.5|) predicts the eigenvalue density at x = 0 to be p(0) = B jn, thereby implying 
B > 0. This condition is satished only for g > g^r = A^^/4. When the coupling g decreases 
below ^'^/4 the one-matrix model undergoes a phase transition and the eigenvalue density 
p{x) develops a two-cut structure, 


with 


7rp(x) = ^^{x^ -b^){a^ -x^) 


+ 2y/g) 
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(2.7) 


( 2 . 8 ) 


At the critical point, where g = gcr = 9^ the solutions for the strong and weak couplings 
coincide. They yield the eigenvalue density which vanishes at x = 0, 


T^Pct{x) = 


g^x‘^ 

2^2 


1 




(2.9) 


For small x this critical density behaves as Pct{x) cx |x|'^ with 5 = 2. The value of the 
critical index 5 is a universal property of the one-matrix model—it does not change even 
when we modify the matrix model potential by any generic higher order terms. In fact, 
5 is related to the so-called string susceptibility 7str = 1 — 5, a quantity which plays an 
important role in string theory. 

Let us see how the critical density, and hence the universal property Pcr{x) oc x^, 
can be extracted from the hydrodynamic picture. The Euler equations which govern the 
motion of the one-dimensional fluid are 













where P is the fluid pressure, p —its density and v —the velocity. For the special case 
when P and p are related by (|^) these equations can be simplifled and even integrated. 
To this end, introduce the complex valued function f — v + inp. In terms of / the system 
( p.lOl) reduces to the complex Hopf equation 


dt dx 


= 0 


( 2 . 11 ) 


which has the following formal solution 


/(a:.*) = /o[a: - 


( 2 . 12 ) 


The function /o(a;) — f{x,t = 0) represents the requisite initial data for the Hopf equation. 
In our problem, one specifies only the real part of /o given by ( p..5|) while the imaginary 
part Im/o(x,0) = 7 rp{x) must be calculated to satisfy the boundary conditions (|1.6|) . 

Due to the time reversal symmetry of Euler equations, imposing ( |1.6|) is equivalent to 
demanding that the velocity at f = 1/2 vanishes identically, v{x,t = 1/2) = 0. Therefore, 
the value of the Hopf function f{x,t) at t = 1/2 ought to be purely imaginary, /(x, 1/2) = 
inpii 2 {x). As a result, the complex equation 

i'Kpii2{x) = /o[a: - \inpii2{x)\ (2.13) 


must have a real-valued solution for Pi/ 2 {x). 

In the limit of large e the critical density (|^ ) and, more generally, (|2.5|) or (p.7|) , 
do give rise to a real pi/ 2 - For example, if p{x) = Pcr{x) the corresponding initial Hopf 
function equals 


fo{z) = v{z) -t- mpcr{z) 

,2 


^ 4 4 3 

- Y" + T" - — 


1 - 




(2.14) 


where we have used (|1.5|) , substituted g = Pcr — and continued /o(^) to z > 2^2/p. 
Now, we shall see in a moment that the typical values of z = x — z7rpi/2(a:)/2 appearing 
in equation (|2.13| ) are of order one, ^ ~ 0{1) S> 1//U. Therefore, we can expand foi^z) for 
large p and finite z. 


/o(^) — —z -h 0 {p 


- 2 ^ 


(2.15) 


which, together with (p.l3|) implies 


'^Pi/2{x) = 2 Vl - -t- 0{p ^). 


(2.16) 














A careful inspection of the above computation shows that equation ( 2.16|) holds for any 
eigenvalue density satisfying the Riemann-Hilbert equation (|2.2|) . Actually, it is simply a 
consequence of the normalization condition f p(x) dx = 1. Indeed, it is easy to compute 
fo(z) given any general U(x), 


fo{z) = ^U'{z) -z + mp{z) 


-z + 


p{y) dy 

z-y ' 


(2.17) 


For an eigenvalue distribution localized around zero in a small interval of width oc l//i this 
formula can be expanded in 1/z to yield ( ^.151) and, consequently, (|2.1ti| ). 

As a result, the shape of 7rpi/2(a:) in the large p approximation says very little about 
the properties of the underlying matrix model. In fact, if we started with Tipi/ 2 {x) = 
2 ^/T — x"^ at t = 1/2 and let it evolve until t = 1 according to the Euler equations, 
that distribution would have collapsed to a delta function. The actual properties of the 
matrix chain are therefore encoded not in the semicircular shape of ( 2.16 ) but rather in 
tiny perturbations over it. These perturbations are strongly amplihed as the distribution 
collapses and produce the nontrivial eigenvalue density p{x, t = 1) = p{x) at the hnal point 
of time evolution. 

Quite remarkably, the effect of small perturbations becomes very transparent if we 
consider the analytic structure of Pi/ 2 (a^)- One discovers that pi /2 has, in general, several 
different regular branches. Equation ( p.l6| ) approximates, of course, only one of them—the 
“physical” branch. However, for e > ecr it happens to be another “hidden” branch that 
actually controls the shape of p{x) = p{x,t = 1). To see this take fi/ 2 {x) = i7rpi/2{x) as 
the initial condition for the Hopf equation and evolve it backwards in time thus obtaining 
/(x, t = 0) — fo{x). The functions fi /2 and /o are related by a formula analogous to (|2.13| ) 


h(x) = /l/2[*+ po(l)]. 


(2.18) 


Now we can substitute x — —iy and use the explicit form of /o(x) from (|2.14|) to derive 
the following parametric representation for 7rpi/2 


7rpi/2(f6 = Fiy) 

^ = -y + ^Fiy) 


where 


F{y) = y + ^y + ^y^-^y^ 




(2.19) 


( 2 . 20 ) 

















Let us emphasize that this representation applies to the special (but most interesting) case 
when the matrix model coupling g has been tuned to the critical value g = gcr{^)- 

The plot of 7rpi/2(*0 ^ function of ^ is shown in hg. 2 and clearly reveals two 

distinct branches. For small ^ one of them behaves as 7rpi/2(*^) ~ 2-^/1 + which is 
simply an anaytic continuation of the physical branch (|2.16|) . The values of the parameter 
y corresponding to this branch are of order one, j/ ~ 1 3> 1///, so that F{y) = y + 1/y + 

Another branch of 7rpi/2 —call it the “hidden” branch—is obtained when y -C 1/^. 
Then for large g we have F{y) k, y{l + IF) and therefore 


7rpi/2(*6 = + 


+ .... 


( 2 . 21 ) 


To demonstrate the importance of this hidden branch let us imagine that we have expanded 
the critical eigenvalue density (|2.9| ) in powers of x. As it turns out, such an expansion is in 
one-to-one correspondence with the expansion of the hidden branch in powers of That 
is to say, given a power series for 7rpi/2(tO (|2.21|) and using equation ( ^.18|) one can 
restore, order by order, the power series for pcr{,x). 



Fig. 2 : The midway density 7rpi/2(*0 lor imaginary arguments x = and the 
graphic solution of equation 7rpi/2(iC) = C/'^- The plot clearly shows two branches 
of 7rpi/2. 


The same would not by any means be true for the semicircular physical branch (|2.16|) . 
There an expansion around zero has nothing to do with how Pcr{x) behaves at small 
X. From this point of view, it is quite satisfactory that the complicated and nonlocal 
perturbations of the physical branch can be summarized in a simple Taylor expansion of 
another branch. 















Of course, to get the particular graph shown in hg. 2 we have used the one-matrix 
model limit e 3> 1. However, we shall see in the next section that the plot of T^Pi/2{iO 
remains qualitatively the same for any e > 1. Therefore, the conclusions derived from such 
a plot—say, the existence of two branches for 7rpi/2(i^) —are not at all restricted to the 
large e limit. 

For a hnite e the hidden branch of 7rpi/2(*^) is no longer approximated by (|2.21| ). 
Instead, it has a certain—yet unknown—Taylor expansion 

7rpi/2(*6 = ai(e)^ + a2(e)^^ ... (2.22) 

all the Taylor coefficients being functions of e. As we advertized above, this expansion can 
be used, together with ( |2.18|) , to reconstruct the potential U{x) and the eigenvalue density 
p{x) = p(a;, t = 0). To this effect we consider the expansion 


/o(x) == fix f 2 x‘^ fsx^ + ... 


(2.23) 


and solve (|2.18|) order by order in powers of x. This yields 


/i = - 

/2 = i 

/s = 


2ai(e) 
ai(e) - 2 
8a2(e) 
[ai(e) -2]^ 
32ai(e) 


(2.24) 


16 a 3 (e) 


[ai(e)-2] [ai(e)-2]' 


and so forth. The real part of fo{x) is related, by virtue of ( |1.5|) and (|2.4|) , to the parameters 
and g. In particular, fi — -|- 1) = —2e^ which hxes ai(e) in terms of e 


ai(e) 


2e^ 


£ 2-1 


(2.25) 


We see that the hidden branch develops a singularity and ceases to exist at e = Ccr = 1. 
This is precisely the location of the Kosterlitz-Thouless phase transition in the inhnite 
matrix chain with the potential ( p.4|) . Remarkably, the same critical value of €cr can be 
derived independently by studying the singularities of the e < 1 phase and by a variety of 
other methods iiil- 

Therefore, the two-branch structure is a characteristic property distinguishing the 
e > 1 phase of the inhnite matrix chain. It contrasts sharply with the one-branch behavior 
found in the e < 1 phase of the same theory [^. Indeed, for e < 1 the small x properties 


















of the critical eigenvalue density can be restored from the physical 7rpi/2(a;) alone, without 
continuing it to any other branches. 

The presence of the second branch has a very interesting reflection in how the droplet 
evolves with time. It results in the formation of a shock—a state where the spatial deriva¬ 
tive of the fluid density becomes inhnite at a certain point. To see how this happens 
consider the time dependence of the droplet density at the origin p{x = 0,t) = d{t)/n. We 
can determine the function d{t) from an implicit relation analogous to (|2.12| ) 

f{x,t) ^ fi/ 2 [x + Tf{x,t)] (2.26) 

T being the shifted time r = 1/2 — t. Setting here a; = 0 and remembering that v{x = 
0 ,t) = 0 at any t, we obtain a transcendental equation on d{t) 


d{t) = 'Kpi/2[iTd{t)\. 


(2.27) 


This equation can be easily solved graphically using the plot of 7rpi/2(*0 shown in hg. 2. 
One simply writes ^ = rdit) so that 


7rpi/2(*6 = ^ (2.28) 

and hnds the intersection point of the plot for 7rpi/2(iO with the straight line of the slope 
1/r. At t = 1/2, when r = 0, this yields a nonzero d given by the value of the physical 
branch at x = 0. As we go backwards in time away from t = 1/2, the fluid density at 
the origin changes. In the specihc case depicted in hg. 2 it hrst grows and then starts 
decreasing until, at a certain moment r = Ter, the density at x = 0 vanishes. This critical 
time Ter is determined by the slope of the hidden branch at the origin 
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Xcr 


d'npi/2{.i0 

d^ 


e=o 


(2.29) 


For any r > Ter or, equivalently, for t < ter = 1/2 — Ter equation (|2.28| ) admits only the 
solution d{t) = 0. Note that, as a consequence of (|2.25| ), the critical time is a function of 
the lattice spacing ter = l/(2e^). 

As a result, the droplet density at the origin remains strictly zero at 0 < t < ter and 
becomes positive in the interval ter < t < 1/2. Of course, the densities at larger times 
t G [1/2,1] are easily obtained using the time rehection symmetry d(t) = d{l — t). The 
shock-type conhguration forms at t = ter and t = 1—ter- This can be understood by solving 
equation ( |2.26|) for small, but nonzero, values of x and times t close to ter- At such t the 
graphical solution to equation (|2.28|) yields a value of ^ very close to zero. Consequently, 












we can use the Taylor expansion of 7rpi/2(*0; given by (|2.22|) , as input for equation (|2.26| ) 
thereby obtaining 


a2(e) 


Tr 


T - Tc 
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p{x, t) cx < 


a? (6) 


2a2(e) 
Ter - r 
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t a2(e)r^rc 


+ a2(e) 


Tc 


Tex - T 


when t < ter 

when t = ter 
B 

x'^ when t > ter 


(2.30) 


for X ^ 0 and \t — tcr| ^ 1- A typical plot of p{x,t) at various times is shown in £g. 3. 
When t < ter the fluid flows smoothly towards x = 0, its velocity increasing in time. 
While t approaches ter, more and more fluid keeps coming in, causing the density gradient 
{d/dx)p{x = 0,tcr) to blow up. This shows as a cusp in the plot of p(x,tcr)- After that, 
the fluid starts accumulating at x = 0, and the velocity x(x,t) becomes smooth again. 



Fig. 3 : The droplet density p{x, t) at various stages of evolution: t < ter (dotted 
line), t = tcT (solid line) and t > ter (dashed line). The plot of p{x,tcr) exhibits a 
“shock” at X = 0. 


Nothing similar ever occurs in the e < 1 phase. There at g = Qeri^) the droplet density 
behaves as p{x, t) oc a{t)\x\ + 0 {x‘^) with a{t) being an inhnitely smooth function of time. 
In that phase p(x = 0, t) = 0 for all t G [0,1]. 

One could say that the quantity ( = pi/2(0) is an “order parameter” for our theory. 
Indeed, C = 0 identically at any e < 1. On the other hand, when e > 1 the values of ( are 


















manifestly positive. Although such an order parameter does not describe the spontaneous 
breaking of any symmetry, it does parametrize a change in the group representation struc¬ 
ture of our model. This shall be demonstrated in section 4. But before that let us use the 
intuition derived from the large e limit to construct an exact solution of the infinite matrix 
chain at e > 1. 


3. Exact Solutions of the Matrix Chain for Lattice Spacings Greater than One 

In this section we shall obtain a class of potentials U{x) for which the chain model can 
be solved exactly. We shall see that the critical index 5 corresponding to these potentials 
equals 5(e) =2 at any e between one and infinity. 

The solutions presented below are not as explicit as exact solutions are usually ex¬ 
pected to be. Neither the interaction potentials nor the eigenvalue densities corresponding 
to these solutions can be written in terms of elementary or special functions. Instead, they 
are determined by complicated transcendental equations, although the critical indices for 
the “transcendental” solutions are still calculable explicitly. 

In principle, to produce such exact solutions is not hard. Choose an arbitrary, but 
properly normalized pii 2 {x). Imagine a droplet which at t = 1/2 is at rest and has this 
density. Evolve this droplet backwards in time up to t = 0 and read off the fluid velocity 
v{x) = v{x,t = 0) together with the density p{x) = p{x,t = 0). The result is an exact 
solution—given by p{x) —to the matrix chain model where the potential U{x) can be 
found from (|1.5| ). 

The most difficult part of this program is to adjust ^ 1 / 2 ( 3 :) so that the potential U (x) 
comes out physically reasonable. Let us require that U'(x) be a smooth and bounded from 
below double-well potential with a hump at a; = 0, 

[/(x) = —^x^ + 0(x'^). (3.1) 

Here, as before, pi^ = 2(2e^ — !)■ 

The potentials in this class are qualitatively similar to the pure quartic potential (|2.4| ). 
Therefore, due to the phenomenon of universality, we can expect that the critical exponents 
of the chain model with these potentials are strictly equal to the critical exponents of the 
pure quartic model. This hypothesis shall be confirmed in section 5 where we investigate 
the pure quartic case using a systematic strong coupling expansion. 

A carelessly chosen pi /2 would most likely yield the potential with < 2, which 
corresponds to e < 1. To generate U{x) with > 2 the analytic structure of pi /2 must be 



very special. This structure can be inferred from the results for the e S> 1 limit obtained 
in section 2. That is, let us choose pi/ 2 {,x) so that 

(a) as an analytic function of x, it has at least two branches (not counting the branches 
different only by a minus sign); 

(b) one of these branches is an even function of x, real valued for real x (this shall be the 
“physical” branch representing the actual droplet density at t = 1/2); 

(c) the “physical” branch is normalized, / pii 2 {x)dx = 1; 

(d) the second branch of pi /2 has the Taylor expansion given by (|2.22|) . 

Although these requirements are quite restrictive, there are many functions obeying 
all of the conditions (a)-(d). As an example, consider the function r[x) = [7rpi/2(a;)]^ 
specihed by the following parametric representation 


r{i^) = R^{ip) cosh(«;(p) 
^ = i?(<p) sinh(K(p) 


(3.2) 


where (/? G [0,7r/2] and R{(p) is an even analytic function of Lp such that R{'k/2) = 0. For 
concreteness, let us choose 

^max 

R{^) = Rn cos"'(^ (3.3) 

n=l 


with arbitrary real coefficients Rn- 

It is easy to check that the plot of r{i^) has two branches as in £g. 3 and 
“physical” branch of r is even. Furthermore, the Taylor expansion of the “hidden 
has the form which agrees precisely with (|2.22D , 


that the 
” branch 


r(fO = 


cosh(7rft;/2) 

sinh^(7rK/2) 


1 + e 


K cosh(7r«;/2) 
.Ri sinh^(7rK/2) 


2 — tanh^ 




(3.4) 


so that condition (d) is also satished. The values of the physical branch at real x can be 
found from (|3.2| ) using the substitution ip — —ip, 




r{x) — cos{kp) 


Rn cosh*^ 


V 


n=l 


^max 

X —sm{K,p) R„cosh"^77. 

n=l 


(3.5) 


As we see, the function r(x) vanishes at ry = ±7 t/2k. Therefore, the support of 7rpi/2(x) = 
\/r{x) is restricted to the finite interval [—aJmax, ^^max] with x^ax = = ±tt/2k), and 

therefore the normalization integral J pi/ 2 {x)dx is hnite. The value of this integral can 
















be adjusted to equal one by rescaling all the coefficients XRn with an appropriately 

chosen A. 

The potential U{x) which corresponds to (|3.2|) must be a nonsingular infinitely dif¬ 
ferentiable function of x. This condition should be enforced at least for those x where 
p{x) 7 ^ 0 —otherwise additional spurious singularities in p{x) could be produced and the 
universality argument (the paragraph after (|3.1|) ) would become inapplicable. 

Technically, U and p are determined by the t = 0 values of the Hopf function f{x,t) = 
v{xR) +i7rp{xR) which, in turn, is constrained by equation ( p.26|) . Consequently, / solves 
a system of transcendental equations obtained from (|3.5|) by the substitution 

{ X ^ X + rfix, t) 

(3.6) 

r ^ -f{xR). 


Generally, this equation has multiple roots. As r = 1/2 — t increases from 0 to 1/2, some 
of these roots might coalesce producing unwanted singularities in U and then also in p. 
Therefore, the parameters Rn in (|3.3|) should be chosen so that this does not happen for 
any r < 1/2. 

The numerical analysis of such transcendental equations indicates (although we do 
not have a rigorous proof) that we do have some choice of parameters for any e > 1. This 
immediately allows us to determine the critical index 5 in the e > 1 phase of the matrix 
chain. Indeed, using the Taylor expansion (|3.4|) we can read off the coefficients 



cosh(7rK/2) 

sinh^(7rK/2) 


cosh(7rK/2) 

3/2 

2 - tanh^ (—^ 

_sinh^(7rK/2) _ 


L V 2 


(3.7) 


From (|2.25|) we now deduce that the parameter is a function of e, 

cosh(7rft;/2) / 2e^ 


sinh^(7rK/2) 


1 


(3.8) 


and that for any lattice spacing a 2 (e) > 0. Consequently, the eigenvalue density in the chain 
model behaves as 7 rp{x) = Im/(x,t = 0) oc / 2 X^, where /2 > 0, yielding the eigenvalue 
index in the e > 1 phase (f(e) = 2. 


By replacing cosh(K(p) and sinh(K(p) in (|3.2|) with more complicated functions it is 


possible to obtain models where at certain discrete values of the lattice spacing e the coef¬ 
ficient a 2 (e) vanishes. At these e such models exhibit higher order multicritical behavior. 





















Of course, the existence of multicritical points in the inhnite matrix chain should not come 
as a surprise. Indeed, exactly the same spectrum of critical indices occurs in one-matrix 
models which describe the e S> 1 limit of the chain theory. 

In the case of the upside-down quartic potential V{M) = m^M‘^/2 + gM'^/A with 
m? > 0 and ^ < 0 the e ^ 1 limit of the matrix chain becomes a c = 0 theory characterized 
by 5 = 3/2. Using the arguments of this and the previous sections it is possible to show that 
there exist regular potentials leading to 5(e) = 3/2 at any e > 1. Again, the midway density 
Pi/ 2 (x) develops a second branch when the lattice spacing e gets bigger than one, while 
only one branch enters the picture for e < 1. However, to see the second branch one should 
continue Pi/ 2 {x) not to imaginary x — but rather to real positive very large x ^ 1/m. 
Since for real x the density Pi/ 2 {x) must always be real, the corresponding “hidden” branch 
can only have half-integer singularities Pi/ 2 ix) cx (x — xq)'^ with 5 = n -|- 1/2, the lowest 
nontrivial value of 5 being S = 3/2. 


4. The Character Expansion and the Order Parameter 


In this section we shall develop the character expansion for the chain model. We shall 
hnd that the SU{N) representation dominating the chain partition f un ction at large N 
develops a “gap” at the Kosterlitz-Thouless phase transition point, and that the gap size 
is a simple function of the order parameter ( = pi/ 2 ( 0 ). 

There are both conceptual and technical reasons why the character expansion is useful 
in the matrix chain problem. Conceptually, it shows that the Kosterlitz-Thouless phase 
transition and the Douglas-Kazakov-type transitions in two-dimensional QCD and the 
dually weighted graph models all occur through essentially the same mechanism. On a 
technical level, there is some evidence, based on the singularity structure of the matrix 
chain with e < 1, that the exact solution of the chain model with a pure quartic potential 


may belong to the class of elliptic functions |]T^. If this is true, the character expansion 
could be just the right computational tool to derive such a solution. 

To construct the character expansion for the inhnite matrix chain, let us go back to 
the partition function (|1.2|) and diagonalize each matrix A4n = UnAnU/ with diagonal 
An = diag(Ai,„ .. .Ajv,n)- Then, express the Hermitian matrix measure dA4n in terms of 
U’s and lambdas, 

dA4n A (Aj^) d\\ n • • • d\j<[^n dUn 


where A (An) is the Van der Monde determinant of eigenvalues Xi^n and dUn refers to the 
Haar measure on the unitary group SU{N). Furthermore, it is useful to introduce the 
matrices 14, = U^_|_iUn. 






Given this notation the partition fnnction of the inhnite matrix chain can be written 
in the form 

Z=j nA'(An)rfAi,„...dA^,„cfG„exp|iVtr ^ - G(A^)] |. (4.1) 


The matrices Vn represent the angnlar degrees of freedom so typical of all mnltimatrix 
models. Fortunately, since the one-dimensional lattice does not have closed loops, all Vn 
are mutually independent and we can easily integrate them out. To do this one simply 
changes variables from {Un} to {Vn} with the result 

Z=j n A"(A^)dAi,^...dA^,„dK.exp|iVtr ^ |. (4.2) 

Note that no Jacobian arises when we pass from U to V. 

In the large N limit the partition function ( [4.2|) can be simplihed even further. Indeed, 
imagine that all integrals over Vn have been done. Then for N ^ oo the remaining 
integrals over the eigenvalues are dominated by an n-independent saddle point Xi^n = \- 
Remarkably, the same saddle point values of A^ occur in a simpler integral 


Zek= / dV A2(A)(iAi .. VAk expjAftryAytA- 


(4.3) 


where A = diag(Ai,..., Ajv)- To prove this assume 


dV exp [iVtr(l/Al/U')] = exp [iV^F(A, A')] 


with a certain F'(A, A'), and write down the saddle point equations for both integrals. 
These equations happen to be the same, and so are their solutions. In fact, the equivalence 
between Z and Zek is just a very special case of the Eguchi-Kawai reduction theorem 

0 . 

There are several ways to integrate out the matrix V in (|4.3|) . One of them—based on 


the large N asymptotics of the famous Itzykson-Zuber integral |]^ — leads eventually 
to the hydrodynamic picture ( p..4|) - (|I.b| ). Another useful method is the character expansion 
||I4|| . Instead of trying to do the integral over V explicitly, we expand the integrand in a 
series of SU{N) characters and then use the formula 


dVxRiVAV^B) = 2 ^xr(A)xr(B). 


( 4 . 4 ) 
















Here Xr{U) stands for the character of a unitary matrix U in the representation R-, 
is the dimension of that representation, while A and B are arbitrary matrices. The hnal 
output of a character expansion is a sum over representations of SU{N) which can usually 


be done by the saddle point method [M |. 

When dealing with large N limits, it is often more convenient to expand in the char¬ 
acters of U{N) rather than SU{N). The U{N) representations are labelled by N integers 
Ri A ■ ■ ■ A the character of a matrix Ai in any such representation being given by the 
Weyl formula 


Xr(M) = 


det II Ap'^ I 

~A(xy 


(4.5) 


As usual, the numbers Ai,..., Aw denote the eigenvalues of the matrix A4, whereas the 
(strictly ordered) integers k equal k = Ui + N — i. 

To compute the integral ( |4.3|) let us expand the term containing W = VAV^A in a 
character series 

= (4.6) 

R 


^NtrW 


The coefficients of this expansion are easily found from the orthogonality of characters. 


Cfl =y dVe 

m 


NtrV^,* 


hi...I n'. 


Xr{V) 

Nl ]X'di +■■■+’'n)- N(N -1)/2 ^ 


(4.7) 


Now we can apply ( [4.4|) to do the unitary group integration in (|4.3|) . Throwing out a 
constant overall factor and using the formula for the dimension of a representation R = 
{h, • • •, In} 

= n ^ (4.8) 


i<j 


this hnally yields 


2ek = 


OO 




h,...,lN=0 


dAi...dAwA2(A) x|^(A)e-'^^''^A). 


(4.9) 


I) are 


At large N both the sum over h, ■ ■ ■ and the integral over the eigenvalues in 
dominated by saddle points. Formally, these saddle points are dehned as the values of 
h,... Jn and Ai,... Aw that maximize the integrand of (|4.9|) . In fact, the saddle point for 
{Ai} is already known to us. It is given by the same eigenvalues that dominate the original 
matrix chain integrand (|1.2|) at N ^ (yo and which are described by the density p{X) in 


(ED 












At the saddle point in the representation space the valnes of li are all of order N. 
Therefore, it shall be convenient to introdnce the qnantities hi = U/N which, in the large 
N limit, condense to form a smooth distribution pi{h). Since all of the /’s are different 
from each other, this density can never be greater than one. 

To maximize the integrand of (|4.9|) with respect to /i,..., /jv we hrst use the Stirling 
formula ^ exp[Nhi{l — log hi)]. Also, the large N asymptotics of the character 

Xi?(A) is given by oc exp{A^^S[p/(h), p(A)]} where the functional S[p/(h), p(A)] has 

a hnite large N limit. In view of this notation, the saddle point equation for the Vs can be 
written in the form 


^ ^ S[pi(h),p(A)] = logh. 


dh Spi{h) 


(4.10) 


The equations needed to compute S are provided by the method of |]^. However, we 


should be particularly careful as some of the Xi (in fact, half of them for an even U{X)) 
have a negative sign. Then the sign of an expression like Ap = Xp^‘‘ oscillates when N 
increases and its large N behavior is not well dehned. To circumvent the difficulty, imagine 
that Ai,...,A 7 v /2 are all positive and that the negative ones are just —Ai,..., —A 7 v/ 2 - 
Furthermore, assume that among h half are odd (call them /°) and half are even (denoted 
/|.) Certainly, if the set of all (odd and even) /’s converges to form the distribution 
Pi{h), the same distribution would also describe the densities of both {/“} and {/|} taken 
separately from each other. 

For such symmetric distributions of {A^} and {h} the Weyl determinant factorizes 
into a product of two {N/2) x {N/2) determinants Q 

(4.11) 


det II Ap" II = (—2)^/^ det ||a!^p|| det ||a!^p| 


where X+p {p = l,...,A^/2 ) are the N/2 positive eigenvalues Ai,..., Ajv/ 2 - Both 
of these smaller determinants exhibit identically the same leading large N behavior 
exp{N^V[pi{h), p{X)]} with A"* = N/2 being the new determinant size. 

To hnd V we hrst represent A+p in the manifestly positive form A+p = exp(6*p/2) so 
that X_/p = exp(A*6*p/ip. The density of thetas a{9) is, of course, related to p(A_|_). Since 
N^a{9) d9 = Ap(A_|_) (iA_|_, one gets 


a 


{9)=e^/^p{e^/^). 


(4.12) 


Now we can write down the equations which hx V. Let us introduce the two functions 


f G+i9) 

[g_(/i) 


d SV 
mda{9) 
d SV 
dh Spiih) 


+ inaij)) 

■ - mpiih). 


(4.13) 








(4.14) 


Then the equation on T) is simply |11][19 


G+[G-{h)] = h. 


This complex-valued functional equation is equivalent to two real-valued constraints which 
determine the two unknown functions {d/d9)ST>/da{9) and {d/dh)ST>/dpi{h). 

On the other hand, the functional derivatives of V can be extracted from the saddle 
point equations. In this way, equation ( [I.10|) yields 


d 5V 
dh Spi{h) 


\ogh. 


(4.15) 


We can also replace Xi —> ±exp(6>j/2) in (|4.9|) , maximize it with respect to 9i and get the 
saddle point equation for the thetas. Remembering that N = 2W, we derive 


d 5V 

m5a{9) 


le^/2c7'(e^/2). 


(4.16) 


Finally, let us combine (|4.13|) , ( |4.15|) and (|4.16|) together with equation ( |4.14|) . The result 
can be conveniently expressed in terms of an auxiliary function H{z) = U'{z)/2 + iTip{z) 
and reads 


H 


y/h 


exp 


inpi{h) \ 

2 J 




(4.17) 


Very remarkably, this equation leads to a very simple formula connecting pi to pi/ 2 - Indeed, 
from (|1.5|) and ( |2.18|) we have the following relation between pi /2 and H, 


i'^Pi/2 


i [H{x) -t- x) 


H[x) — X. 


(4.18) 


Now let us substitute x — v^exp[—i7rp;(h)/2] and use ( [4.17|) . This produces, at last, the 
formula 


7rpi/2 


Ai/h 


cos 


/ T^pi{h) \ 

V 2 )\ 


= 2 Vh sin 


2 ; 


(4.19) 


which is our hnal result. 

Equation (|4.19|) is essentially kinematical. Indeed, it does not involve the matrix model 
interaction potential in any way0. However, it does indicate that the descriptions of the 
inhnite chain model based on the hydrodynamic picture and on the character expansion 
are completely equivalent. 


^ Equation ( 4.19 ) can also be derived in an “explicitly kinematic” fashion following the ideas 
of Douglas 











































Fig. 4: The Young tableau densities pi{h) for e < 1 (left) and for e > 1 (right). 
The presence of a gap in the Young tableau for e > 1 shows as a plateau on the 
graph of pi{h). 


To be rigorous, equation ( [4.19|) holds only for those h where pi{h) is strictly less than 
one. Indeed, if the Young tableau density reaches its upper bound pi{h) = 1 the saddle 
point equation (|4.15|) ceases to apply and the whole derivation becomes invalid. Practically, 
this allows us to determine whether there is, in fact, an interval of h where pi{h) = 1. If 
yes, then the size of this interval h* equals the “gap width” in the corresponding Young 
diagram as shown in hg. 1. 

Since h* is just on the border of the region where (|4.19|) applies, we can use that 
formula with h = /i* and pi{h^) = 1. This yields ( = 7rpi/2(0) = 2\/h^ thereby proving 
equation (pT7|). In other words, the gap size, which appears to be a natural “order parame¬ 
ter” for the phase transition in the character expansion language is related to C, —an order 
parameter in the hydrodynamic picture. 

Therefore, in the critical regime h* = 0 for e < 1 while /i* = h*(e) > 0 for any 
e > 1. This picture of the Kosterlitz-Thouless phase transition is completely identical 
to the mechanism of the Douglas-Kazakov transition in large N two-dimensional QCD 
BH. Indeed, both in the matrix chain model and in QCD 2 the density pi{h) develops a 
plateau for large lattice spacings (respectively, strong coupling) while pi{h) < 1 everywhere 
for small lattice spacings (or weak couplings.) And, to complete the correspondence, both 
transitions are induced by the topologically nontrivial states—vortices and instantons— 
that dominate in the strongly coupled phases of these two theories. 



















5. The Chain Model with a Pure Quartic Potential 


In this section we shall constrnct and explore a systematic strong coupling expansion 
for the inhnite matrix chain with a quartic interaction potential. The expansion we shall 
consider is in powers of 1/e^ or, equivalently, 

For large fj.^ the successive terms in the character expansion ([4.91) become smaller 
and smaller. Indeed, any character Xi?(A) is a homogeneous polynomial in the eigenvalues 
Xi. The degree of homogeneity for such a polynomial equals the total number of boxes 
in the corresponding Young tableau, njj — ni + ... + un = li + ... + In — N{N — l)/2. 
Consequently, if for large fi we approximate the eigenvalue density p(A) by the one-matrix 
model result (|2.5|) , the typical width of such distribution will be Atyp ~ l//r and the value 
of the character xh(A) will be suppressed by a factor Ij. In other words, to expand in 
1/^ up to the order 1 /we should choose only those representations in ([4.9|) which have 
not more than ur boxes. The characters in such representations shall be polynomials in 
trA, trA^,..., trA"^^ and would effectively generate corrections to the one-matrix potential 
U (A). That is to say, they will perturb the right hand side of the Riemann-Hilbert equation 
( p.2| ) by the terms of order at most 

A one-matrix model with such modihed potential has, in general, the solutioni 


Tip{x) = Q{x)\ 1 




(5.1) 


where Q{x) is a polynomial of degree ur — 2 and —the “renormalized” distribution 
width [^. Obviously, if the original potential U{X) was even, so shall be Q{x). The 
coefficients of Q could, in principle, be determined by computing the characters, evaluating 
the corrections to U{X) and solving the resulting one-matrix model. 

Fortunately, there is a more economical way to organize this computation. Let us use 
the fact (following essentially from (|2.12|) ) that for any inhnite chain model the 

functions 

G±{z) = (5-2) 


(with p{z) being the exact eigenvalue density) obey the functional equation 


(5-3) 

Very remarkably, this equation allows one to compute Q{z) much easier and quicker. 

^ Since we are interested in the critical regime, the coupling constant g is always adjusted to 
equal 5cr(e). For the uncorrected one-matrix model this would imply g = 








The left hand side of equation (|5.3|) contains the eigenvalue density at a complex 
point, p[Q-[z)]. Consequently, it is important to choose the correct analytic branch of p in 
both Qj^ and Q-. Since it is more convenient to work with real variables, we shall enforce 
equation (|5.3|) in the region of large positive z > 2^2 !where i'Kp{z) becomes real. The 
correct sign choice for p{z) in that region yields 

g^(z) = g^(z) s g(z) = (5.4) 

The reason is, p{z) has a square root branch cut between Xmax = 2\/2lp^: and —Xmax = 
—2\f2lp^. If a real 2 were positioned on the upper edge of this cut the value of 0 -(z) 
would have a negative imaginary part. In evaluating p[Q_{z)] this would take us below 
the lower edge. But the values of p{z) on the lower edge are opposite to the values of p 
on the upper edge. When z is continued to z > 2^1 p^ the branch of p^z) entering 
should be continued from above the cut while p{z) in Qj^ should be continued from below 
the cut. This compensates for the sign difference in front of p{z) in (|5.2|) and leads to the 
result (1^). 

At inhnitely large p, which corresponds to the one-matrix model case, we get 

e+(z) = e-(2) = I + o(Tj (5.5) 

so that equation (|5.3|) is trivially satished. The 1 /corrections can now be derived by 
perturbing this solution. Since we are interested in the critical regime where p{x) cx 
the polynomial Q{z) must vanish at z — 0, and we can seek a solution of the form 


Giz) 


-2; H—2 

2 2 




(5.6) 


The coupling g = gcr{p) and the coefficients of P{z) should be computed from equation 
( b.3|) together with the normalization requirement J p{x) dx = 1. 

Technically, it is easier to expand in powers of 1/^^ rather than l/p'^. This is perfectly 
legitimate given that, to the leading order in our large p expansion, p = p^.. In other words, 
we shall seek corrections to the one-matrix model relations p = p^ and pcr = p^j^-, 



(5.7) 
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(5.8) 


where we cut the expansion at order 1/. For the polynomial P{z) we assume 
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Note that, to 0{l/n‘^^), the degree of P{z) equals 2(M — 2), fully consistent with (|0|). 

To determine the coefficients bn-, Qn, ei and as,i we impose Q[Q{z)] — z + 0{ii~‘^^) and 
require the normalization of the density to be 1 + Practically, the normalization 

requirement amounts to asking that the coefficient of I/ 2 : in the large 2 ; expansion of (|5.6| ) 
equal one. The resulting system of equations can be solved recursively using the following 
two stage procedure. First, one expresses 6 m+i, fi'M +2 and ag^i with s + z = M — 1 in terms 
of Cm and the coefficients of lower order. Then cm is determined by the normalization 
equation, but only at order M+2. In this manner, all the coefficients in (|5.6|) are determined 
unambiguously. 

The l/fi expansion described above can be easily carried out to very high orders. Let 
us present here just the hrst few terms of this expansion. 
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(5.9) 
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The information about the critical properties of the theory is encoded in the value of 
P{z — 0). Indeed, the critical eigenvalue density p{z) behaves at small z as np{z) ~ a{p^)z‘^ 
with 


9 


a(/) = ^P{z = 0) 
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1 + 


1 ^ 

■“ 9? 


(5.10) 


If a{p^) is hnite and nonzero, the critical index 5 equals 5 = 2. On the other hand, a 
vanishing ai^p?) would indicate a higher order multicritical point. Also, one could encounter 
a singularity where a{p^) goes to inhnity. This would mean that for small x the density 
p{x) vanishes slower than and the asymptotics p{x) cx does not apply. In fact, we do 
expect such a singularity to occur at the Kosterlitz-Thouless point e = 1 or, equivalently, 
p? = 2. At that point [|^ 


p(x) cx 


X 


log[l/(A|x|)] 
which is indeed more singular than p(x) cx x^. 


(5.11) 



















Fig. 5 : The plots of a(//^) obtained from the strong coupling expansion with 
M = 8,10,12 and 14 (counting from lower left to upper right.) It is seen that 
a(/r^) does not vanish and exhibits a singularity around = 2. 


Therefore, the coefficient a(/r^) should blow up at = 2. In fig. 5 we plot a(/r^) 
computed from the strong coupling expansion up to These plots show that 

certainly does not vanish. Therefore, no higher order multicriticalities occur for e > 1. Fur¬ 
thermore, does blow up in the region around = 2 . The Fade approximants reveal 

only one singularity in that region which we can therefore associate with the Kosterlitz- 
Thouless phase transition. Of course, to determine the critical value of numerically is 
very hard, as the strong coupling expansion becomes unreliable when gets less than 
~ 3. However, when combined with the theoretical results found in section 3, the nu¬ 
merical data become quite conclusive. Indeed, the universality arguments suggest that 
the difference between the “transcendental” potentials of section 3 and the pure quartic 
potential could only lead to a higher order multicritical point in the quartic case. But such 
points are reliably excluded by the data from the strong coupling expansion. At the same 
time, the exact position of the Kosterlitz-Thouless point can be derived from equation 
(p.25|) independently of any numerical simulations. 









To summarize, the strong coupling expansion, taken together with the exact results, 
indicates that the inhnite chain model with a pure quartic potential exhibits the critical 
index 5 = 2 for any e > 1. 


6. Conclusions 

In this paper we have explored the e > 1 regime of the inhnite random matrix chain. 
We have calculated the critical exponents and found that, for e > 1, they are the same 
as the critical exponents in one-matrix models. As a consequence, the critical behavior of 
the strongly coupled two-dimensional 0(2) model on a random planar lattice turns out to 
coincide with the critical behavior of pure two-dimensional gravity. 

From the technical viewpoint, we found that most properties of the chain model are 
encoded in the midway density Pi/ 2 {x) which enters the picture through the hydrodynamic 
representation. We have seen that the changes in the branch structure of pi /2 govern the 
Kosterlitz-Thouless phase transition and explain the change in critical exponents. 

Finally, we constructed the character expansion for the matrix chain. Rather remark¬ 
ably, in the character expansion language the mechanisms of the Kosterlitz-Thouless phase 
transition and the Douglas-Kazakov-type transitions are precisely the same. We derived 
an exact relation between the “Young tableau density” of the representation dominating 
the character sum at large N and the hydrodynamic density pi/ 2 - 

There are, however, a number of problems that remained open. First of all, it would 
be desirable to engineer an exact solution reproducing the critical behavior of the matrix 
chain on both sides of e = 1. According to (|5.11J ), the requisite midway density pi /2 would 
have to contain a logarithmic singularity and, for e > 1, have the two-branch analytic 
structure. We are conhdent that such a solution exists and can be constructed with our 
methods. 

Furthermore, it would be very interesting to interpret our results in terms of another 
theory describing the two-dimensional 0(2) model—the matrix quantum mechanics on a 
circle. There, too, the representation structure changes at the Kosterlitz-Thouless point, 
although in a somewhat different way. Since the chain model and the compactihed matrix 
quantum mechanics are dual to each other, there might exist a correspondence between 
the large N limits of their character expansions. 

In this respect, the matrix theory on a circle resembles the instanton description of 
QCD 2 . Indeed, two-dimensional QCD on a sphere has two dual descriptions. One of 
them is the character expansion given by the Migdal-Rusakov formula The other 

represents the partition function of QCD 2 as a sum over all instanton sectors, and can be 
obtained from the character expansion by Poisson resummation [pT|. 




At large N the representation dominating in the character expansion formula has 
0{N) boxes and, depending on the phase, may or may not have a gap. The Young 
tableaus corresponding to such representations are shown in £g. 1. On the other hand, 
for N ^ (yo the instanton expansion is also dominated by a single sector. The N charges 
parametrizing the dominant instanton all strictly equal zero at weak coupling but become 
nonzero integers of order N at strong coupling |]^ [|^ . 

This is very similar to what happens in the matrix quantum mechanics on a circle. 
There the vortexless phase is dominated by the trivial (singlet) representation, with zero 
boxes in each Young tableau row [|5| [^ . In the other phase, where the vortices are signif¬ 
icant, the dominant representation has 0{N) boxes. We see that the row lengths of the 
dominant Young tableau in the matrix quantum mechanics behave the same way as the 
instanton charges in QCD 2 . 

All of this suggests that, by analogy with QCD 2 , a simple Poisson resummation of 
the character expansion formula ( |4.9|) could tell us a lot about compactihed matrix quan¬ 
tum mechanics. If true, such analysis could be extremely useful, as the matrix quantum 
mechanics on a circle so far appeared intractable by any known computational method. 
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